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Abstract 

In this paper, we employ the bifurcation theory of planar dynamical systems to 
investigate the exact travelling wave solutions of a generalized Degasperis-Procesi 
equation ut — Uxxt + ^uUx + 'yiu — Uxx)x = ^UxUxx + uUxxx- The implicit expression of 
smooth soliton solutions is given. The explicit expressions of peaked soliton solutions 
and periodic cuspon solutions are also obtained. Further, we show the relationship 
among the smooth soliton solutions, the peaked solitons solution and the periodic 
cuspon solutions. The physical relevance of the found solutions and the reason why 
these solutions can exist in this equation are also given. 
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1 Introduction 

Recently, Degasperis and Procesi [1] derived a nonlinear dispersive equation 

Ut - Uxxt + 4:UUx = SUxUxx + UUxxx (1-1) 



which is called the Degasperis-Procesi equation. Here u{t, x) represents the 
fluid velocity at time t in the x direction in appropriate nondimensional units 
(or, equivalently the height of the water's free surface above a flat bottom). 
The nonlinear convection term uUx in Eq.(l.l) causes the steepening of wave 
form, whereas the nonlinear dispersion effect term Su^Uxx+uUxxx = {\u^)xxx in 
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Eq.(l.l) makes the wave form spread. Eq.(l.l) can be regarded as a model for 
nonlinear shallow water dynamics [2]. Degasperis, Holm and Hone [2] showed 
that the Eq.(l.l) is integrable by deriving a Lax pair and a bi-Hamiltonian 
structure for it. Yin proved local well-posedness to Eq.(l.l) with initial data 
uq G H'^ (M) , s > 3 on the line [3] and on the circle [4] . The global existence of 
strong solutions and weak solutions to Eq.(l.l) are investigated in [4]-[10]. The 
solution to Cauchy problem of Eq.(l.l) can also blow up in finite time when 
the initial data satisfies certain sign condition [7]- [10]. Vakhnenko and Parkes 
[11] obtained periodic and solitary-wave solutions of Eq.(l.l). Matsuno [12,13] 
obtained multisoliton, cusp and loop soliton solutions of Eq.(l.l). Lundmark 
and Szmigielski [14] investigated multi-peakon solutions of Eq.(l.l). Lenells 
[15] classified all weak traveling wave solutions. The shock wave solutions of 
Eq.(l.l) are investigated in [16,17]. 

Yu and Tian [18] investigated the following generalized Degasperis-Procesi 
equation 

— Uxxt + 4:UUx = SUxUxx + UUxxx ~ JUxxx, (1-2) 

where 7 is a real constant, and the term Uxxx denotes the linear dispersive 
effect. They obtained peaked soliton solutions and period cuspon solutions 
of Eq.(1.2). Unfortunately, they didn't obtain smooth soliton solutions of 
Eq.(1.2). 

In this paper, we are interesting in the following generalized Degasperis- 
Procesi equation 

Ut - Uxxt + 4:UUx + 7(ti - Uxx)x = ^UxUxx + UUxxx, (1-3) 

where 7 is a real constant, the term Ux denotes the dissipative effect and the 
term Uxxx represents the linear dispersive effect. Employing the bifurcation 
theory of planar dynamical systems, we obtain the analytic expressions of 
smooth solitons, peaked solitons and period cuspons of Eq.(1.3). Our work 
covers and supplements the results obtained in [18]. 

The remainder of the paper is organized as follows. In Section 2, using the 
travelling wave transformation, we transform Eq.(1.3) into the planar dynam- 
ical system (2.3) and then discuss bifurcations of phase portraits of system 
(2.3). In Section 3, we obtain the implicit expression of smooth solitons and 
the explicit expressions of peaked solitons and periodic cuspon solutions. At 
the same time, we show that the limits of smooth solitons and periodic cusp 
waves are peaked solitons. In Section 4, we discuss the physical relevance of the 
found solutions and give the reason why these solutions can exist in Eq.(1.3). 
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2 Bifurcations of phase portraits of system (2.3) 



We look for travelling wave solutions of Eq.(1.3) in the form of u{x,t) = 
ip{x — ct) = ^p{^), where c is the wave speed and ^ = x — ct. Substituting 
u = ip{^) into Eq.(1.3), we obtain 

- ap' + ap'" + Aifip' - Sip'cp" - ifiif'" + -fip' - -fip'" = 0. (2.1) 

Integrating (2.1) once we have 
ip"iip - c + 7) = - (c - 7)v. + 2(/.2 _ (^/)2^ ^2.2) 

where g is the integral constant. 

Let y = Lp', then we get the following planar dynamical system: 



dy ^ g-{c--f)ip + 2tp'^ -y'^ 



(2.3) 



with a first integral 

H{^, y) = {^-c + 7)' {y' -^^-g)=h, (2.4) 



where h is a constant. 

Note that (2.3) has a singular line = c — 7. To avoid the line temporar- 
ily we make transformation = ((/? — c + 7)(iC- Under this transformation, 
Eq.(2.3) becomes 

^ = ((/?- c + 7)?/ 

(2.5) 



dy 



g - {c - 'y)(p + 2(p^ - y 
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System (2.3) and system (2.5) have the same first integral as (2.4). Conse- 
quently, system (2.5) has the same topological phase portraits as system (2.3) 
except for the straight line = c — 7. Obviously, = c — 7 is an invariant 
straight-line solution for system (2.5). 

For a fixed h, (2.4) determines a set of invariant curves of system (2.5). As 
h is varied, (2.4) determines different families of orbits of system (2.5) having 
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different dynamical behaviors. Let M{(fe,ye) be the coefficient matrix of the 
hnearized system of (2.5) at the equihbrium point {ipe,ye), then 



M{iPe,ye 



ye ipe- C + J 



(2.6) 



and at this equihbrium point, we have 

J{^e, ye) = det M((^e, Ve) = "2?/^' " (^e - C + -f){4ipe " C + 7), (2.7) 
p{(Pe, ye) = trace(M((/7e, ye)) = -ye- (2.8) 

By the quahtative theory of differential equations (see [19]), for an equilibrium 
point of a planar dynamical system, if J < 0, then this equilibrium point is a 
saddle point; it is a center point if J > and p = 0; if J = and the Poincare 
index of the equilibrium point is 0, then it is a cusp. 

By using the first integral value and properties of equilibrium points, we 
obtain the bifurcation curves as follows: 

9i{c) = (2.9) 

^2(c) = -(c-7)'. (2.10) 

Obviously, the two curves have no intersection point and (72 (c) < < gi{c) for 
arbitrary constants c 7^ 7. 

Using bifurcation method of vector fields (e.g., [19]), we have the following 
result which describes the locations and properties of the singular points of 
system (2.5). 

Theorem 2.1 For given any constant wave speed c ^ 0, let 



c - 7 ± J(c - 7)2 - 8fi^ 
= for g<gi{c), (2.11) 



y^ = ±^{c-j)^ + g for g>g2{c). (2.12) 
When c = 7, 

(1) if g < 0, then system (2.5) has two equilibrium points (—-^—1,0) and 
-|,0), which are saddle points. 
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(2) if g = 0, then system (2.5) has only one equilibrium point (0,0), which 
is a cusp. 

(3) if g > 0, then system (2.5) has two equilibrium points (0,— y^) and 
(0, ^Jg), which are saddle points. 

When c 7^ 7, 

(1) if g < g2{c), then system (2.5) has two equilibrium points {ipQ,0) and 
{ipQ,0). They are saddle points. 

(1) ifo-i, then < -\{c--i) < i(c - 7) < c - 7 < (/j^. 
(a) ifc< 7, then v?o < c - 7 < i(c - 7) < -^(0-7) < (^0 ■ 

(2) if g = g2{c), then system (2.5) has three equilibrium points ((/?o,0), 
(v9q , 0) and (c — 7, 0). (c — 7, 0) is a cusp. 

(i) ifc>-f, then v^o = - 7) < |(c - 7) < c - 7 = (v?o , 0) zs a 
saddle point, while {(Pq,0) is a degenerate center point. 

(li) if c <-i, then (/Jq = c - 7 < i(c - 7) < -|(c - 7) =ipQ. (v?o > 0) is a 
degenerate center point, while ((/Jo',0) is a saddle point. 

(3) if g2{c) < g < gi{c), then system (2.5) has four equilibrium points 
i(Po,0), (v5o,0), (0-7,2/0") and {c-'y,y^). (c-7,?/^) and {c--f,y^) are two 
saddle points. 

(i) if a > ^, then (^q < |(c — 7) < < c — 7. {(f^,0) is a saddle point, 
while {ipQ,0) is a center point. 

(a) if c < J, then c — j < ip^ < |(c — 7) < c/Jq . ((/7q , 0) is a center point, 
while ((/?(]', 0) is a saddle point. 

Specially, when g = 0, 

(i)ifc>'y, then the three saddle points {ipQ,0), (c — 7,2/(7) and (c — 7,2/0") 
form a triangular orbit which encloses the center point {ipQ,0). 

(a) if c < 7, then the three saddle points {ipQ,0), (0 — 7,2/5") and {c — j,yQ) 
form a triangular orbit which encloses the center point {(pQ,0). 

(4) if 9 = 9i{c), then system (2.5) has three equilibrium points (^^,0), 
(c— 7,2/0") (ind {c—j,yQ) . (^^,0) is a degenerate center point, while (c— 7,2/0") 
and (c — 7, 2/0") are two saddle points. 

(5) if g > gi{c), then system (2.5) has two equilibrium points (c — 7,2/0^) 
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and (c — 7, i/q), which are saddle points. 



Corresponding to the case c = 7 and the case c 7^ 7, we show the phase 
portraits of system (2.5) in Fig.l and Fig.2, respectively. 




(a) (b) (c) 

Fig. 1. The phase portraits of system (2.5) (c = 7). (a) (7 < 0; (b) g = 0; {c) g > 0. 

3 Solitons, peakons and periodic cusp wave solutions 

Theorem 3.1 Given arbitrary constant c 7^ 7, let ^ = x — at, then 
(1) when < g < gi{c), 

(1) if c> ^, then Eq.(1.3) has the following smooth hump-like soliton solu- 
tions 

/3i((/;+)=/5i(<^)e-l«l for ^^<^<ipt, (3.1) 

(a) if c < ^, then Eq.(1.3) has the following smooth valley-like soliton 
solutions 

/32(<^) = /32(^-)e-l«l for i^^<^<^t, (3.2) 

(2) when g = 0, Eq.(1.3) has the following peaked soliton solutions 
(^=(c-7)e-l«l, (3.3) 

(3) when g2{c) < g <0, 

(i) if '-f, then Eg. (1.3) has the following periodic cusp wave solutions 
u{x, t) = (fsix -at- 2nT) for {2n - 1)T < x - at < (2n + 1)T, (3.4) 
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(a) if c < 'y, then Eq.(1.3) has the following periodic cusp wave solutions 
u{x, t) = (f^ix -ct- 2nT) for {2n - 1)T < x - at < (2n + l)r, (3.5 



where 



n , . ^ (2 y/y?^ + + /- + + /i) (yp - ^-^Y^ 

^^^^^ ^ (2v/y^ + miLp + ma + 2y? + mx)[}f> - ^tY^ 

+ niiip + 1712 + + rusY'^ 



, 3(c-7) + \/(c-7P-8g 

'1 2 ' ' 



3(c - 7)^ - % + 5(c - 7)v/(c - 7)^ - 8^ 
8 

[c - 7)2 -Ag + 3(c - 7)\/(c-7)2-8^ 



(3.9) 



h = 2 , ^3.10) 

3(c - 7) - v/(c - 7)' - 8^7) ^^^^^ 

mi = , (3.11 

2 ' ^ ' 



3(c - 7)' - 4^? - 5(c - 7) v/(c - 7)^ -^9 
m2 = r , (3.12) 



ms = , (3.13) 



«i = > (3-14) 
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-(c- 


-7)- 


v/(c-7)2- 


-8(?, 







02 = ^ , (3.15) 

f3.16) 



h2 = -{c--i) + d{c--iY-Sg, (3.17) 



3(c-7) + v/(c-7)'-8<7 
CKi = , (3.18) 



2^(0 - 7)2 - 8(? + 3(c - 7) V^- 7)' - ^9 
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Ot2 



-3(c-7) + v'(c-7)'-8^ 



2 J(c - 7)2 - 8(7 - 3(c - 7)v/(c-7)'-8(? 



¥^3(0 = + /_el«l for < (/^3 < c - 7, 

(/74(0 = + ^-e"'^' /or c - 7 < < -v^, 



± 



T 



c-7± \/(c-7)^ + ^ 



-2^) - ln(2L 



(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 



= -(c - 7) ± ^ - - % T 2 V - 7)' T (c - 7)V(c - 7)' - 8f7, (3.24) 



ip^ and are as in (2.11). 

Before proving this theorem, we take a set of data and employ Maple to dis- 
play the graphs of smooth solion, peaked soliton and periodic cuspon solutions 
of Eq.(1.3), see Fig.3-Fig.7. 

Proof. Usually, a solion solution of Eq.(1.3) corresponds to a homoclinic orbit 
of system (2.5) and a periodic travelling wave solution of Eq.(1.3) corresponds 
to a periodic orbit of system (2.5). The graphs of homoclinic orbit, periodic 
orbit of system (2.5) and their limit curves are shown in Fig.8. 

(1) When Q < g < gi{c),c > 7, system (2.5) has a homoclinic orbit (see 
Fig.8 (a)). This homoclinic orbit can be expressed as 



c — 'J — ip 



for Pq <p < Iff. 



(3.25) 



Substituting (3.25) into the first equation of system (2.3) and integrating along 
this homoclinic orbit, we obtain (3.1). 

When < 51 < (7i(c), c < 7, we can obtain (3.2) in similar way. 

(2) When (7 = 0, c > 7, system (2.5) has a homoclinic orbit that consists of 
the following three line segments (see Fig. 8(c)). 

y = ±f for <ip< ift: (3-26) 
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and 

= c - 7 for - yjic-'yy + g <y< \/ {c - -f)"^ + g . 



(3.27) 



Substituting (3.26) into the first equation of system (2.3) and integrating along 
this orbit, we obtain (3.3). 

When 5f = 0, c < 7, we can also obtain (3.3). 

(3) When (72(c) < g < 0, c > ^, system (2.5) has a periodic orbit (see Fig. 
8(e)). This periodic orbit can be expressed as 



y = ±yf^ + 9 for v^<V5<c-7, (3.28) 

and 

(p = c--f for - \J{c-^Y + g <y< ^J{c--iY + g. (3.29) 

Substituting (3.28) into the first equation of system (2.3) and integrating along 
this periodic orbit, we obtain (3.4). 

When g2{c) < (7 < 0, c < 7, we can obtain (3.5). ■ 

Remark 3.1 From the above discussion, we can see that when g < 0, g ^ 0, 
the period of the periodic cusp wave solution becomes bigger and bigger, and the 
periodic cuspon solutions (3.4) and (3.5) tend to the peaked soliton solutions 
(3.3) . When g > 0, g ^ 0, the smooth hump-like soliton solutions (3.1) and 
the smooth valley-like soliton solutions (3.2) lose their smoothness and tend 
to the peaked soliton solutions (3.3). 



4 Discussion 

In this paper, we obtain the solitons, peakons, and periodic cuspons of 
a generalized Degasperis-Procesi equation (1.3). These solitons denote the 
nonlinear localized waves on the shallow water's free surface that retain their 
individuality under interaction and eventually travel with their original shapes 
and speeds. The balance between the nonlinear steepening and dispersion 
effect under Eq.(1.3) gives rise to these solitons. 

The peakon travels with speed equal to its peak amplitude. This solution is 
nonanalytic, having a jump in derivative at its peak. Peakons are true solitons 
that interact via elastic collisions under Eq.(1.3). We claim that the existence 
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of a singular straight line for the planar dynamical system (2.3) is the original 
reason why the travelling waves lose their smoothness. 

Also, the periodic cuspon solution is nonanalytic, having a jump in deriva- 
tive at its each cusp. 
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Fig. 3. Smooth hump-like soliton solutions of Eq.(1.3)( c = 2, 7 = 1). (a) g = 0.12 
(b)g=0.1; (c) g=0.075; (d) g=0.05. 
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Fig. 4. Smooth valley-like soliton solutions of Eq.(1.3).(c = 1, 7 = 2). (a) g 
(b) g=0.1; (c) g=0.075; (d) g=0.05. 
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(a) (b) 
Fig. 5. Peaked soliton solutions of Eq.(1.3). (a) c = 2, 7 = 1; (b) c = 1, 7 = 2. 
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(d) (e) (f) 

Fig. 8. The orbits of system (2.5) connecting with the saddle points. 
(a)0 <g < gi{c),c > 7; (b) < 5 < gi{c),c < 7; (c) 5 = 0,c> 7; (d) g = 0,c< 7; 
(e) 52(c) < g < 0, 7; (f) ff2(c) < 5 < 0, c < 7. 
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